Intermittent outgassing through a non-Newtonian fluid 
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We report an experimental study of the intermittent dynamics of a gas flowing through a column of 
a non-Newtonian fluid. In a given range of the imposed constant flow rate, the system spontaneously 
alternates between two regimes: bubbles emitted at the bottom either rise independently one from 
the other or merge to create a winding flue which then connects the bottom air-entrance to the 
free surface. The observations are reminiscent of the spontaneous changes in the degassing regime 
observed on volcanoes and suggest that, in the nature, such a phenomenon is likely to be governed 
by the non-Newtonian properties of the magma. We focus on the statistical distribution of the 
lifespans of the bubbling and flue regimes in the intermittent steady-state. The bubbling regime 
exhibits a characteristic time whereas, interestingly, the flue lifespan displays a decaying power- 
law distribution. The associated exponent, which is significantly smaller than the value 1.5 often 
reported experimentally and predicted in some standard intermittency scenarios, depends on the 
fluid properties and can be interpreted as the ratio of two characteristic times of the system. 
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PACS numbers: 05.45.-a, 05.40.-a, 47.57-s, 83.60.Rs, 91.40.Hw 



I. INTRODUCTION 

A large number of physical systems driven far from 
equilibrium exhibit an intermittent behavior, where quiet 
phases (off-state) randomly alternate with periods of 
strong activity (on-state). Experimental examples in 
physical systems include electronic devices [l[ , spin- wave 
instabilities |2i] , plasma produced by a gas discharge [1] , 
liquid crystals [J, |5|, and nanoscopic systems such as 
single molecules nanocrystal quantum-dots 0], or 
nanoparticles diffusing through a laser focus Q. Inter- 
mittency has also been reported in other disciplinary 
fields, like in the example of the sleep-wake transitions 
during sleep @ . A standard characterization of the in- 
termittency phenomenon is the distribution of time spent 
in the ojf-state. Interestingly, in most of the experimen- 
tal studies reported so-far [l|, H S i, H, H, 0, Hi, this 
distribution exhibits a power-law regime with an expo- 
nent close to 1.5, with a cut-off at large times. This value 
of the exponent is well-understood theoretically through 
the so-called scenarios of on-off intermittency p^[Tll.[l3| 
and type-Ill intermittency [ij, [3l , which both predict a 
power-law distribution with an exponent 3/2. 

A Newtonian fluid flowing through a non-Newtonian 
material is alike to present different escaping regimes 
[m as well as an intermittent activity, as previously ob- 
served both in physical [l^, [l^] and geophysical [H, [3| 
contexts. In the laboratory, experiments have been con- 
ducted in the case of air emission through an immersed 
granular layer . These studies evidenced different air- 
outgassing regimes through the non-linear medium at a 
constant flow rate. In particular, it has been observed 
that gas emitted at the bottom of an immersed granu- 
lar layer can cross the whole system through ternpqrary 
channels, which later spontaneously collapse [H, [l3|- In 
a geophysical context, volcanologists reported observa- 
tions of the degassing process which interestingly exhibits 
similar temporal patterns [l8|. In this case, the gas. 



released by magma during its rise toward the free sur- 
face, plays the role of the Newtonian fluid and magma 
the role of the non- Newtonian material [l^, [2l|: inside 
the volcano conduit, dissolved gas gives birth to bubbles, 
which ascend and burst intermittently at top of the vol- 
cano vent [l9|. What remains puzzling is the existence 
of spontaneous changes between different activities {e.g. 
discrete bursts or continuous puffing [12]) for which the 
non-Newtonian pro perties of magma could play an im- 
portant role [ii,!!!,!!!]. 

At the laboratory scale, two kinds of experiments have 
been previously conducted in order to model the volcanic 
degassing- activity. On the one hand, experiments were 
designed to account for the potential effects of the geome- 
try of an underlying magma-chamber . These studies 
were limited to the case of Newtonian fluids and, be- 
sides, did not consider the issue of the time spent in each 
of the regimes. On the other hand, recent laboratory 
experiments performed on real rhyolitic melts close to 
their natural conditions have revealed that during their 
rise, bubbles deform and coalesce under shear [2^ lead- 
ing to channel- like bubble networks [l^]. Those struc- 
tures were proposed as a key ingredient in the magma 
degassing process at depth [23|, but not directly related 
to the non-Newtonian rheology of magma. 

In this paper, we report an experimental study of the 
intermittent dynamics of a gas rising through a column 
of a non-Newtonian fluid. We unravel that this out- 
gassing process can be considered as a two-states pro- 
cess : either it occurs through bubbles rising indepen- 
dently (the bubbling regime, on-state), or through a gas 
channel which crosses tortuously the whole system {the 
winding-flue regime, off-state). Let us emphasize here 
that such an intermittent process, defined in terms of a 
two-states process (bubbling or channel- like regime), is 
different from the intermittent regime reported in purely 
Newtonian fluids, as described in [2^. Indeed, a New- 
tonian solution cannot sustain the formation of a flue 
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(in contrast to non- Newtonian fluids [23|) and, in this 
case, intermittency can be related to the formation of 
antibubbles [1^. Here, we report statistical data con- 
cerning the irregular oscillations between the bubbling 
and the winding flue regimes. Specifically, we measure 
the statistics of the time spent in both the bubbling and 
flue regimes. We observe that the lifespan of the flue 
is characterized by a power-law distribution with an ex- 
ponent significantly smaller than the often reported ex- 
ponent 3/2. This non-standard exponent turns out to 
depend slightly on the rheological properties of the fluid. 
Using a simple stochastic two-states model, we relate this 
exponent to the ratio of two characteristic times of the 
system, that we estimate from the experimental data. 

These results might be of interest to the geophysics 
community as the magma, charged with bubbles, exhibits 
a flow-threshold and a shear-thinning behavior [l^ [l^, 
[2ll . |30| . Indeed, non-Newtonian rheological properties 
are likely to explain the existence of different outgassing 
processes (including bubble networks) and to govern the 
associated nontrivial statistics. 



II. EXPERIMENTAL PROCEDURES 
A. Experimental setup 

The experimental setup (Fig. [T|) consists of a vertical 
plexiglass tube (inner diameter 74 mm, height 270 mm) 
filled with a non-Newtonian solution. Air is injected at 
the bottom of the fluid column through an air tank, par- 
tially filled with water. The system thus makes it possi- 
ble to tune the inner free volume V of the tank (ranging 
from 60 to 800 cm^) by changing the inner amount of 
water and to inject humid air. Besides, above the fluid 
column, a small container filled with water maintains a 
saturated humidity level at the free surface. This way, 
one avoids any significant drying of the sample during the 
experimental time. A mass-flow controller (Bronkhorst, 
Mass-Stream series D-5111) introduces air inside the tank 
at a constant flow rate Q (from Qmin = 0.17 cm^/s 
to Qmax = 1-74 cm'^/s). Air is then injected at the 
bottom of the fluid column through a nozzle (diameter 
d = 2.0 mm). Rigid tubes (typical diameter, 8.0 mm) 
insure the connection between the tank and the column. 
We measure the variation of the overpressure SP inside 
the chamber with the help of a low-pressure differential- 
sensor (Honeywell S&C, 176PC28HD2) connected to a 
multimeter (Keithley, 196). A C-I--I- program records the 
pressure data (to within 5 Pa) through the GPIB inter- 
face over long time-durations (typically a few days) at 
5 Hz. As we will see, monitoring the pressure in the tank 
chamber makes it possible to probe the degassing activ- 
ity while direct visual observation is obstructed as the gel 
solution progressively flUs up with bubbles (Sec. IIIC")) . 

The fluid consists in a diluted solution (10% in mass of 
distilled water, unless otherwise specifled) of a commer- 
cial hair-dressing gel {Gel coifFant, fixation extra forte, 




FIG. 1; Sketch of the experimental setup 

mass-flow controller 



Air, from a 

is injected, at a constant flow rate 
Q, in a chamber (volume V) connected to the bottom of the 
fluid column (height h) . By means of a presure sensor P, we 
measure the overpressure 5P inside the chamber as a function 
of time. In the picture on the right-hand-side, one can notice 
that the density of bubbles trapped into the gel strongly de- 
pends on the altitude (The scale is given by the inner diameter 
of the plexiglass tube, 74 mm.) 



Auchan). This choice is mainly justified by the fact 
that one can easily be supplied with large quantities of 
fluid, reproducible mixtures are rather easy to prepare 
f3l]| and they are stable in time if well prevented from 
drying. Some of the non-Newtonian rheological prop- 
erties of the fluid are reported in the appendix A. In 
particular, the solution, when free of bubbles, is shear- 
thinning and we estimated the viscoelastic characteristic 
time Tc — l/7c — 1700 s, compatible with the ratio of 
the viscosity at low shear-rate 770 — 6.10^ Pa.s to the 
estimated flow threshold CTc — 35 Pa. 



B. Preliminary observations 

The gel solution is poured inside the tube and a con- 
stant flow is imposed. Qualitatively, one observes that 
the solution starts filling up with small-sized bubbles. 
These small bubbles, which remain trapped because of 
the fluid threshold, are further advected by the large bub- 
bles subsequently emitted at the nozzle and progressively 
a vertical gradient pointing toward the top of the column 
takes place (Fig. [T]). After this transient regime, whose 
duration depends on the flow rate Q, the fluid column 
reaches a constant height h, significantly larger than the 
initial height Hq. We observe that there exists a fiow rate 
Qf, which can be experimentally defined as the most effi- 
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cient flow rate to fill the solution with trapped bubbles. 
The transient which lasts a few hours for Q ~ can 
last days to weeks if Q significantly differs from Q^. Spe- 
cial attention has been paid to the characterization of the 
bubble gradient and of the associated fluid rhcology. De- 
tails are reported in the appendix A. Finally, note that 
all the data reported in this article were obtained after 
the bubble gradient was established. 



C. Experimental protocol 

We pointed out that small bubbles remain trapped in 
the fluid, which might lead to an evolution of the rheo- 
logical properties through time (aging). Thus, in order 
to ensure the reproducibility of the experiment, we chose 
the following experimental protocol : the gel solution, 
initially free of bubbles, is poured inside the tube; the 
flow rate is set to the smallest accessible value, Qmin', 
the overpressure SP is then recorded for, at least, 3 days. 
Subsequently, the flow rate is increased by steps of am- 
plitude SQ up to the maximum value, Qmax- For each 
value of the flow rate Q, the pressure P is again recorded 
for, at least, 3 days. Then, the same process is appHed 
to decrease the flow rate Q down to Qmin- The whole 
process will be referred to as a flow rate cycle. For all the 
experiments reported here, we chose 6Q ~ Qmin, corre- 
sponding to about ten data points over the experimental 
range of Q. 



III. EXPERIMENTAL RESULTS 
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FIG. 2: Typical pressure variation SP{t) vs time t - 
At a constant flow rate, the typical signal from the pressure 
sensor exhibits a spontaneous alternation between rest and 
activity periods. Rapid drops of the overpressure mark the 
emission of successive bubbles at the bottom of the column 
whereas an almost-constant overpressure corresponds to a flue 
connecting the injection hole to the free surface. The analysis 
consists of measuring the lifespan ti, (resp. r/) of the bubble 
(resp. flue) regime. The figures from top to bottom display 
the signal at different time scales {ho = 70 mm, V = 800 cm'^ 
and Q = 0.35 cm^/s). 



A. Overall behavior of the system 

We now turn to the study of the pressure recorded in- 
side the tank-chamber. We report SP, the variations of 
the pressure around its mean value, high-pass filtered at 
1 mHz. One immediately notices that the pressure ex- 
hibits two distinct behaviors (Fig. [2]), which emphasizes 
rest and activity periods. 

Thus, a first crucial result is the following: for a given 
flow rate, the system spontaneously alternates, in an ir- 
regular manner, between two states. Either cusped-like 
bubbles are emitted at the bottom of the gel column, 
rise through it and burst at the top independently from 
one another (the bubbling regime); the corresponding 
pressure variation 6P exhibits successive rises and drops, 
the latter corresponding to the formation of the bubbles. 
Or, a channel, previously reported as a 'bubble sausage' 
[29| , connects the bottom air-entrance to the column top 
(winding Que regime); the corresponding SP is almost 
constant. 

The formation of the flue is linked to the non- 
Newtonian properties of the hair-gel solution and, as we 
shall see, to the thixotropy. A rising bubble produces 
a shear which locally leads to a decrease of the fluid 
resistance to flow, thus creating an easy path through 



the fluid column (S^- This process takes place if the 
time-difference between two successive bubbles is shorter 
than the typical relaxation time scale of the fluid. Once 
formed, the flue, subjected to a creep flow under the ef- 
fect of gravity, slowly deforms and eventually collapses, 
leading the system back to the bubbling regime. These 
two mechanisms explain why the system switches quickly 
between two very different states. 

As a first global result, we report the fraction of time, 
$b, spent by the system in the bubbling state at a given 
flow rate (Fig. [3]). Below a certain threshold Q*, the 
system experiences the bubbling-regime, whereas above 
Q* the system naturally alternates between the bubbling 
and winding-flue regimes (We note that Q* ~ Qb, which 
is discussed in the appendix B). We point out that the 
two regions, above and below the transition, are not sym- 
metric: indeed, for Q < Q* , we only observe the bubbling 
regime (the system is unable to form a flue) whereas, in 
the case Q > Q*, even for the largest flow rate, the flue 
always eventually tweaks and collapses (finite lifespan). 
The critical flow rate Q* is independent of the column 
height liQ and of the chamber volume V but depends on 
the rheological properties of the fluid. For instance, the 
minimum flow rate necessary to sustain an open flue in- 
creases when the fluid threshold is decreased by dilution 
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FIG. 3: Fraction <J>f, vs. imposed flow rate Q - The 

fraction of time spent in the bubbling regime, $5, clearly 
exhibits a transition at a given flow rate Q* which dras- 
ticaUy depends on the fluid rheology. Fitting curves, sig- 
moids, are given as eye leads. The gray area marks the range, 
Q G]Q*jQ^]i in which the duration of the bubbling regime 
can be deflned (associated with •, only) [(symbol: height ho 
mm, volume V cm"^). Added-water content 10 %: (•: 70, 530); 
(o: 130, 530); (8: 70, 800). Added-water content 15 %: 
(□: 70, 530)]. 



and we assessed Q* ~ 0.3 cni-^/s for 10% of added water 
and Q* ~ 1 cm^/s for 15% (Fig.©. Finally, note that, 
even if the channel can always collapse at large flow rate, 
one can identify a flow rate Q"*" ~ 1.5 cm'^/s above which 
the emission of a single bubble leads to the formation 
of a new channel (flow rate Q above the gray region in 

Fig. ID. 

From now on, because of the binary behavior of 
gel column at a given flow rate Q, we shall account 
for the complex dynamics of the system by consid- 
ering the overpressure signal as a temporal sequence 

r (1) (1) (2) (2) (n) r ;i • -O 

{Tj: ,Tj: ',...,Tj: ' , '\ of thc succcssivc jiue- 

(resp. bubbling-) regime lifespans, tj (resp. r^). We 
aim at estimating the probability distributions Pb{Tb) and 
Pf{Tf). To do so, we display estimates of the correspond- 
ing cumulative probability distributions Fi,{t\j) (Fig. [3]) 
and Ff{Tf ) (Fig.[6l), defined as 

pAr)dT (1) 

where x stands for 6 or /. We sort the sequence {ri*''} 
into {tx''^} such that 

ri(i) > ri(^) > . . . > (2) 
and plot j/n versus t^*'"''' [s^ . 




FIG. 4: Cumulative probability distribution Fi,{ti,) - 
The time Tb spent in the bubbling regime exhibits an expo- 
nential distribution with a characteristic time (Fig. 
We report data obtained when increasing and decreasing 
the flow rate Q [(symbol, Q cm^/s). Top: (x, 0.347); 
(o, 0.521); (■, 0.695); (A, 0.869); (T, 1.04); Middle: (0, 1.22); 
(♦, 0.782); (o, 0.608); (M, 0.434). Bottom: (x, 0.347); 
(o, 0.521); (■, 0.695); (▼, 1.04). (Added-water content 10%, 
V = 530 cm^ ho = 70 mm)]. 



B. Bubbling regime 

From the temporal sequence {t^^"''}, we deduce that 
there exists a characteristic time, , to build a flue from 
the bubbling regime. Indeed, the cumulative probability 
distribution Fb(rb) decreases almost exponentially for in- 
creasing values of Tb (Fig. [4]). Thus, the lifespan Tb of the 
bubbling regime presents an almost-exponential probabil- 
ity distribution. We report the associated characteristic 
time as a function of the flow rate Q for different 
values of the chamber volume V (Fig. [5]). First, r^* de- 
creases drastically when the flow rate is increased above 
Q* (Note that is not defined for Q < Q* as the sys- 
tem is constantly in the bubbling regime and that, for 
Q > Q"*", the measurements are difficult to perform be- 
cause of the small number of channel collapses during 
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Q (cm /s) 



corresponds to the upper limit in figure [S] In order to de- 
fine this upper hmit properly, one must consider instead 
the flow rate for t^* to equal the typical time of a bub- 
ble emission. We estimate, in this case, that the channel 
forms immediately, after the emission of a single bubble, 
for Q > 1.5 cm^/s. 

Thus, even if this rough model must be considered with 
caution, the proposed mechanism accounts for the exis- 
tence of a characteristic lifespan r^* of the bubble regime 
in a finite range ]Q*,Q'^] of the imposed flow rate, 
decreasing when the flow rate Q is increased. 



FIG. 5: Time vs. flow rate Q - Values of obtained for 
both increasing and decreasing flow rates are reported. The 
closed symbols denote the values deduced from the data re- 
ported in the figured The open symbols correspond to data 
obtained when the fluid column is stirred after each change 
in the flow rate (□, = 530 cm^). The continuous line cor- 
responds to the equation (|3]) with ro ~ 65 s, Qc = 0.35 cm^/s 
and Q* ~ 0.29 cvo? /s. The horizontal dashed line denotes 
the characteristic time of emission of a single bubble, about 
1 s. The grey region corresponds to Q £ ]Q* ,Q^] [(Sym- 



C. Winding flue regime 



From the temporal sequence {Tj""*}, 



^ ^ , , we observe that, 
for Q G the cumulative distribution Ffijf) of 

the winding-flue lifespan Tf presents a power-law behav- 



bol, V cm^) : (•, 530); (♦, 800). (Added-water content 10%, ior, Ff{Tf) oc " for r/ £ [ri,T2] with ti - 10 s and 
ho = 70 mm)]. T2 ~ 1000 s (Fig. [6]) and a rapid decay beyond the cut-off 

time T2. We report in the flgure [7] the experimental val- 
ues of a obtained for flow rates in the gray region defined 
in the figure [31 We estimate that a — 0.18 ± 0.02, this 
latter value being independent of Q and V . However, a 
depends on the spatial distribution of the small bubbles 
which alter the rheological properties of the fluid (Ap- 
pendix A). This conclusion is drawn by using a slightly 
different protocol: in an experimental series, we stirred 
the solution after each increment in the flow rate and 
estimated an even smaller value a ~ 0.1 in this latter 
case (Fig. [71). We thus observe that a is always signifi- 
cantly lower than the usual value a = O.SLoften reported 
in experimental observations [J [2, d, 0, 0, d, 0, [1| and 
also expected from some of the theoretical intermittency 
scenarios, like the on- Ojff intermittency [lOj HHj ■ Note 
that such a small value of the exponent {a < 1) imposes 
a cut-off at long times, T2, which is necessary for the aver- 
age (r/) to exist. Experimentally, the data correspond to 
an acquisition time of 4 days, which is significantly larger 
than T2, so that we expect the system to have reached a 
stationary state [s^l and the cut-off not to be due to 
the finite experimental time. In some sense, T2 corre- 
sponds to the largest memory time scale during which 
the stresses in the fluid at rest remain correlated, which 
is compatible with the viscoelastic time of the fluid, that 
is to say T2 ^ Tc, thus of the order of 10'^ s (Appendix 
A). 

Hence, the flue lifespan does not exhibit any charac- 
teristic timescale in a temporal window between ti and 
T2 in which a power-law distribution with a small expo- 
nent is observed. These observations, which lie at the 
core of the present article, are thoroughly discussed in 
the section ITVl 



the experimental time.) Second, we observe that is 
almost independent of V. In contrast, depends on the 
distribution of the small bubbles in the fluid. Indeed, 
increases slightly when the ffuid column is stirred each 
time the flow rate is set to a new value (Fig. which 
reveals the sensitivity of the phenomenon to the rheology 
of the fluid. 

The existence of a characteristic time t^* associated 
with the lifespan of the bubbling regime is quite intuitive 
but the mechanism leading to the formation of the chan- 
nel is complex. A full description of the phenomenon, 
which would deserve an extensive study, is beyond the 
scope of the present article. However, we propose a sim- 
ple heuristic argument (Appendix B) which suggests to 
describe the experimental data by the following fitting 
function (Fig. [5l): 



-To In 




(3) 



The key ingredients are a decrease of the flow thresh- 
old (yield stress) due to the shear imposed by the ris- 
ing bubbles and the thixotropy of the fluid charged with 
bubbles. Interestingly, the rheological measurements re- 
ported in the appendix A point out that the small bub- 
bles enhance the thixotropy and thus play here an im- 
portant role. From the interpolation of the experimental 
data, we obtain tq — 65 s, Qc — 0.35 cm^/s and Q* ~ 
0.29 cm^/s (added-water content 10%, V = 530 cm^ and 
V ~ 800 cm'^, ho = 70 mm). Interestingly, requiring 
> in the equation we obtain a maximum fiow 
rate = QcQ* /{Qc — Q*) — 1-7 cm^/s, which roughly 
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FIG. 6: Cumulative probability distribution Ff{Tf) - 
The cumulative probability distribution Ffijf) of the time 
spent in the flue regime exhibits a power-law behavior, up 
to a large-time cutoff T2. Curves have been vertically shifted 
for clarity, and the slope 1/2 characterizing standard on-off 
intermittency is shown for comparison. Note that the data 
correspond to an acquisition time of 4 days, which is sig- 
nificantly larger than T2, so that we expect the system to 
have reached a stationary state [sj and the cut-off not to 
be due to the finite experimental time [(Symbol, Q cm'^/s). 
Top: (o, 0.521); (■, 0.695); (A, 0.869); (T, 1.04). Mid- 
dle: (o,1.22); (♦, 0.782); (o, 0.608); (M, 0.434). Bottom: 
(x, 0.347); (o, 0.521); (■, 0.695); (T, 1.04). (Added-water 
content 10%, V = 530 cm^, ho = 70 mm)]. 



IV. DISCUSSION 

As explained above, our main result is the experimen- 
tal observation that the probability density function of 
flue lifespans displays a power-law behavior over a broad 
time window, with a non-standard exponent 1-1- a ~ 1.2, 
which slightly depends on the rheological properties of 
the fluid. This result is of particular interest, since in 
standard intermittency scenarios, the exponent takes a 
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FIG. 7: Exponent a vs. flow rate Q - The values of a ob- 
tained for increasing and decreasing fiow rates (respectively 
closed and open symbols) are almost the same and signifi- 
cantly smaller than 1/2. Stirred samples (crosses) lead to even 
smaller values of a (added-water content 10%, V = 530 cm^ 
and V = 800 cm^, hg = 70 mm). 



universal value, independent of the details of the dynam- 
ics. For instance, in the on-off intermittency scenario, 
the exponent 1 -I- a = 3/2 is directly related to the first 
return time probability of a random walk (or a "chaotic 
walk" in the deterministic case) [l2| . In purely determin- 
istic scenarios, like type-I, type-II and type-Ill intermit- 
tency scenarios [l^ [ij] , the exponents are related to the 
classification of bifurcations in dynamical system theory 
[13, [H, US • The corresponding values of the exponent 
1 -I- a in these scenarios (1/2 for type-I intermittency, 2 
for type-II and 3/2 for type-Ill [llfll) all differ signifi- 
cantly from our experimental values. 

In the present experiment, due to both its non- 
standard value and its observed dependence on the ex- 
perimental conditions, a is likely to be related to some 
rheological or physical properties of the system, rather 
than to some generic theoretical property. It is thus in- 
teresting to relate the measured value of the exponent 
a to some properties of the gel. A full analysis of this 
issue is well beyond the scope of the present article, but 
in order to get some first insights, we propose to describe 
the data with a simple stochastic model with two states, 
namely bubble or flue. In such a framework, we write 
that, after having spent a time in the bubble regime, 
the system has a probability per unit time (or transi- 
tion rate) ri,{Tif} to switch to the flue regime. Similarly, 
we assume that after a time Tf in the flue regime, the 
transition rate to the bubble regime is rf{Tf). In the bub- 
ble regime, the exponential distribution, _F),(Tf,), with the 
characteristic time r^* means that the transition rate rf, is 
constant, independent of t^, and equal to rf, = 1/t^. In 
contrast, in the winding-flue regime, the observed power- 
law distribution Ffij-f) implies that the transition rate 
rf(Tf) is not a constant, and thus that a more complex 
dynamics is at play, even in this phenomenological de- 
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scription. From the relation 



dTf 



-rf{Tf)Ff(Tf), 



(4) 



a power-law behavior Ff{Tf) (x Tj:" is found to corre- 
spond to rf(Tf) = a/r/. Still, the power-law behavior is 
observed only in a window ti < Tf < T2, which suggests 
that, as a first approximation, r/ could be taken as a con- 
stant outside this window. A simple parameterization of 
rf{Tf) leading to such a behavior is given by 



^0(1 + t//ti) 



T2 



(5) 



provided a = ti/tq and ri T2. The resulting distribu- 
tion 



Ff{Tf) 



-Tf/T2 



(6) 



behaves as a power law for t larger than ri, with an 
exponential cut-off at T2. In the opposite limit r/ ^ ri, 
Ff(Tf) can be expanded as 



(7) 



This linear decay of Ff{Tf) at short time can thus be 
used to evaluate the time scale tq. In addition, the time 
scale Ti can be estimated from the crossover between the 
linear and power-law behavior. From the experimental 
data, we assess tq ^ 30-40 s and n ~ 5-8 s. These 
values lead to a ratio ti/tq in the range 0.12-0.27, com- 
patible with the experimental value of a. We can thus 
consistently interpret the exponent a as the ratio of two 
characteristic time-scales: the time ti beyond which the 
winding flue starts to age (roughly the time for the flue 
to reach its stationnary shape) and the time tq typically 
needed for the flue to collapse in the early stationary 
regime. Finally, note that the time tq is of the same or- 
der of magnitude as the time needed for the yield stress 
of the fluid to relax to its asymptotic value Uc after the 
shear induced by the rise of a single bubble (see appendix 
11). 



V. CONCLUSION AND OUTLOOKS 

We have shown that the flow of a Newtonian fluid into 
a non-Newtonian one can exhibit an intermittent activ- 
ity, namely between a bubbling and a winding flue regime. 
The distribution of the ofl-state lifespan displays a power- 
law behavior with an exponent about 1.2, which signif- 
icantly differs from the standard exponent 3/2. A sim- 
ple stochastic two-states model suggests that this non- 
standard exponent could be directly related to the ratio 
of two timescales characterizing the experimental system. 

Our results also contain relevant conclusions in a geo- 
physical context, even if we do not claim at any direct in- 
terpretation. We propose that the volcanic intermittent 



activity observed on the field, which has been hitherto 
explained in terms of constriction of the magma cham- 
ber [2^ or by c hang es in gas flux and/or magma flow 
rate at depth '35] , could be (partly) the results of the 
non-Newtonian properties of magma. Indeed, our simple 
laboratory experiment demonstrates that, even if a con- 
stant gas flow is imposed, one can observe an intermittent 
activity on top of the fluid column. 

A more detailed understanding of this non-standard 
intermittency phenomenon would certainly be desirable, 
both from practical and fundamental viewpoints. An ex- 
tensive study of the dependence of the exponent a on 
the physical properties of the fluid (rheological proper- 
ties, surface tension [3§|, . . . ) is in progress. Moreover, 
we are also currently investigating a two-dimensional ex- 
periment, similar to the present three-dimensional one, 
in which one can monitor the evolution and the collapse 
of the flue. We aim at relating the overall dynamics of 
the system to the local dynamics. 

Finally, we proposed that the minimum ingredients 
necessary for a channel to form and stabilize, and thus 
for the intermittency to be observed, are a shear-thinning 
and thixotropic fluid. However, a detailed study of the 
mechanisms at stake in the flue formation and collapse 
remains to be performed. 
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APPENDIX A: FLUID CHARACTERIZATION 

1. Rheology of the fluid 

The rheological properties of the solution, free of bub- 
bles, have been measured using a rheometer (Bohlin In- 
struments, C-VOR 150) equipped with parallel plates 
(PP-60, gap e = 750 and 1000 /im). Sand paper was 
glued to the plates in order to prevent any sliding of the 
fluid at the walls; At 25°C, we performed both steady- 
shear and steady-stress experiments with an equilibration 
time At ranging from 1 to 25 s (Fig. [5]). The solution is 
shear-thinning and we estimated the associated charac- 
teristic time Tc = l/7c — 1700 s, compatible with the 
ratio of the viscosity at low shear-rate 770 ~ 6 x 10^ Pa.s 
to the estimated flow threshold ctc — 35 Pa (Fig. [5]). 



2. Bubbles size and density 

During the experiments, the small bubbles that remain 
trapped in the fluid column, due to the yield stress, al- 
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FIG. 8: Fluid viscosity rj vs shear-rate 7 - The fluid is 
sliear-tliinning. Tlie vertical gray line indicates the charac- 
teristic shear-rate 7c. Inset: Fluid viscosity rj versus shear 
stress a - The measurements clearly point out the existence 
of a flow threshold CTc — 35 Pa [Open symbols correspond 
to stress-controlled experiments whereas solid symbols corre- 
spond to shear-controlled experiments. Gap e = 750 /im : 
(□, At = 25 s); {<}, At = 1 s); (•, At = 10 s)]. 



ter the rheological properties of the fluid. We report here 
measurements of the size and density of the small bub- 
bles along the fluid column when the system has reached 
a steady state. 

During the first cycle of Q, the state of the column 
evolves as it fills up with small bubbles which get trapped. 
However, after having switched Q twice to Q* during the 
increase and the decrease of the flow rate, the aspect 
of the column remains qualitatively the same, and one 
can assume that the vertical bubble gradient has reached 
a steady state (We give several indications of such an 
assertion in Sect. IIIII ) In order to characterize the verti- 
cal bubble gradient which takes place along the column 
height in the steady-state regime, we perform the fol- 
lowing experiment: from a 15 cm high column which 
has experienced 5 flow rate cycles, we carefully extract 
7 slices, 2 cm thick, corresponding to different altitude z; 
For each slice, we measure the size of a thousand bubbles 
using a macroscope {Leitz, Laborlux) and assess the fluid 
density. The corresponding rheological properties shall 
be reported in the section [X3] 

Let us denote as Pnir) the probability distribution of 
the bubble size r in the slice n. In each slice, the cumu- 
lative probability distribution of the bubble size r, de- 
fined as Fn{r) = f^pn{r')dr' is exponential (Fig. [H]). 
Thus, the bubble size, in each slice, is exponentially 
distributed. The corresponding characteristic radius, 
r* = (0.58 ± 0.04) mm, is almost constant over the col- 
umn height h (Fig. [TUb ) . In order to assess the associated 
bubble concentration, we measured the density of mate- 
rial in each slice (Fig. [TOb). The fiuid density decreases 
roughly linearly from the density of the fluid free of bub- 
bles at the bottom to p ~ 750 kg/m^ at the free surface, 
corresponding to a maximum gas fraction of about 25%. 
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FIG. 9: Bubbles size and density - Left: Cumulative prob- 
ability distribution Fn of the bubble size r in the slice n (The 
grayscales are uniformely distributed from light gray at the 
bottom to black at the top.) The distribution Fn is exponen- 
tial (linear flt, dotted black-line.) We denote r* the associated 
characteristic radius. Right: The pictures, from bottom to 
top, correspond respectively to the slices, z G [0.5, 2.5] cm, 
z € [8.5, 10.5] cm and z £ [12.5, 14.5] cm, respectively. 



The measurements confirm the visual observation of a 
bubble gradient inside the fluid. We estimate that the 
bubble concentation, which equals at the bottom, is 
about 300 bubbles per cm'^ at the top. 



• (a) 



{0} 



-it- 



14 
12 
10 
8 
6 
4 
2 


0.0 0.4 0.8 800 1000 4000 5000 25 35 45 



r„* (mm) 



P (kg/m ) (Pa.s) (Pa) 



FIG. 10: Fluid column properties - We report (a) the 
mean bubble size r*, (b) the density p, (c) the viscosity rjo 
in the limit of low shear-rate 7, and (d) the fluid-threshold 
(Tc along the column height z. Each point corresponds to an 
average over one slice (thickness 2 cm, horizontal dashedlines). 



3. Rheological properties in presence of bubbles 

Using the same experimental protocol as in the section 
lA 11 we determined, for each slice, the rheological proper- 
ties of the fluid, charged with bubbles. We observe that 
the bubbles lead to the appearence of a significant mem- 
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FIG. 11: Viscosity rj vs shear stress a - The rheological 
properties depend on the ahitude z because of the variation 
of the density of the bubbles trapped in the fluid along the 
column height. Moreover, we observe that the data obtained 
for increasing and decreasing shear rate (arrows) do not su- 
perpose: the presence of the bubbles leads to a significant 
memory effect. Inset: Fluid viscosity rj vs shear-rate 7 - The 
measurements clearly point out the memory effect (The sym- 
bols, from light gray to black, correspond respectively to the 
slices, z € [0.5, 2.5] cm, z € [3.5, 5.5] cm, z € [8.5, 10.5] cm 
and z G [12.5, 14.5] cm, where z denotes the altitude. Imposed 
shear-rate, equihbration time At = 15 s). 



ory effect (Fig. [TT|) : the shear stress does not reach the 
steady-state value instantaneously but after a significant 
characteristic time, tq, which we estimate from the ex- 
perimental data to be about 1 niin for the largest bubble 
density (at the top of the column). In contrast, we ob- 
serve that the shear-thinning properties of the fluid are 
not qualitatively altered (Fig. [11] inset). Note finally that 
the yield stress Cc — 25 Pa is roughly constant along the 
column (Fig. [TUB), and smaller that the estimated value 
for the free bubble gel solution (ctc — 35 Pa). 



APPENDIX B: ESTIMATE OF < 

In order to account, at least at a rough estimate, for 
the observation of a characteristic time of the bubbling 
regime, we propose that a channel can form and remain 
stable if the shear-stress imposed by the inner gas flow 
at the channel walls exactly balances the yield stress. 
Initially, the yield stress is large and the gas emission at 
the nozzle necessitates a large overpressure which leads 
to the formation of a bubble once the fluid starts flowing. 
The rise of successive bubbles leads to a decrease in the 
yield stress which allows the emission of gas at a small 
overpressure and, thus, the formation of a thin channel 
in the fluid. Moreover, once the channel has formed, 
there is no flow of the fluid at the walls and the yield 
stress increases back to its maximum value, which tends 
to stabilize the structure. 

On the one hand, between two successive bubbles, the 



fluid is at rest and the yield stress, CTy, relaxes toward 
its asymptotic value ac over the characteristic time tq 
according to 



To CT y 



(Bl) 



where tq is an intrinsic property of the fluid. On the 
other hand, due to the large shear stress imposed by the 
rise of the bubble, ay tends to vanish according to 



0. 



(B2) 



where is here not an intrinsic property of the fluid. In- 
deed, Tb characterizes the destructuration of the fluid un- 
der stress and, for instance, might depend on the rising- 
bubble size. Hence, Tb < tq jSTj . 

Considering the typical time, r = Vb/Q, between two 
successive bubbles {Vb denotes the volume of the bub- 
ble) and the duration of the shear by a rising bubble, t;,, 
according to equations (|B2p and (|Bip . one can estimate 
that, after the rise of n successive bubbles: 



CTc 

T Tb 



1 — exp — n — 

To 



(B3) 



where we assumed that tb t ^ Tb and only kept the 
exponential dependence on n. 

When the yield stress has reached a small enough 
value, the gas emission can occur in a gentle manner, a 
thin gas finger growing from the nozzle. We can estimate 
that the growth of the finger will lead to a stable channel 
provided that the shear stress applied to the wall by the 
inner gas fiow does not exceed the yield stress. Assuming 
a Poiseuille flow inside, one can estimate the shear stress 
at the wall ap = ArjaQ / {Trr^) where rja = 1.6 x 10~^ Pa.s 
denotes the dynamic viscosity of air and Tc the radius of 
the channel. In this framework, the condition ct; 
leads to: 



p 



-To In 



tb To 



(B4) 



where we deflned = ht. This latter equation can be 
rewritten : 



-To In 




where, by deflnition. 

To tb 



and Q* — 



Qc 



l + Qc/Q^ 



(B5) 



(B6) 



with Q 



+ — 



Note that t^ (> 0) is thus predicted to 



4rja ■ "'""^ 'b 

be defined in a finite range Q € ]Q* , Q"*"] of the flow rate 
and, in agreement with the experimental observations, 
not to depend on the height of the fluid column. 

From the interpolation of the experimental data 
(Fig. Qc - 0.35 cmVs, Q* ~ 0.29 crn^/s and 
To ~ 65 s. First, we observe that To is found to be in rea- 
sonable agreement with the estimation obtained in the 
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rheology experiment. Second, from the time difference 
between successive bubbles in the pressure signal (Fig. [2]), 
we obtain the typical bubble-volume Vb — 2 cm^ and, ac- 
cordingly, the typical radius ~ 0.8 cm. We estimate 
the rise time of the bubble along the whole column height 
(8 cm) to be about 1 s and, thus, tb ~ 0.2 s (about 2rb/vb, 
the diameter divided by the rise velocity, Vb — ^ cm/s). 
With this latter value of tb and the experimental Qc, we 
obtain Tb ~ 2.3 s, which is a reasonnable value if one con- 
siders that the shear rate imposed by the rising bubble is 
large, of about 1/tb — b s^^ (Note, in addition, that the 
imposed stress is of about pgrb ~ 80 Pa where p denotes 
the density of the fluid and g the acceleration due to the 
gravity.) From Q+ ~ 1.7 cm^/s and CTc ~ 25 Pa (Fig.fTO]). 
one can estimate the channel radius rc ~ 0.1 mm, which 
is also in agreement with the experimental observation 
of a very thin channel, a fraction of millimeter in diam- 
eter. At last, it would be interesting to understand the 
mechanism that selects the channel radius (competition 
between the elasticity and the surface tension?) but the 
problem is out of the scope of this appendix. 



From the proposed mechanism, one can also under- 
stand why the small bubbles quickly invade the system 
when Q is about Q* (i.e. Qb — Q* , Sec. IIIB|) . Indeed, 
on the one hand, in the bubbling regime, the pinch-off 
of the bubble emitted at the nozzle 38], the coalescence 
of two bubbles in the bulk ([39|, Sec. nT|) and the burst- 
ing of a bubble at the free surface [3l| are likely to leave 
a small satellite (bubble) [i^. We can thus guess that 
each rising bubble leaves one or two small bubbles. On 
the other hand, in the flue regime, the number of small 
bubbles left by the channel collapse is about the column 
height divided by the wavelength that develops (Rayleigh 
instability ^43|). In our experimental conditions, we esti- 
mate from direct visual observation that the wavelength 
is about 1 cm (much larger than the channel diameter) 
and, thus, that each collapse of a channel produces about 
10 small bubbles. The channel collapse is thus much more 
efficient in the production of small bubbles. By definition 
of Q* , at the transition, the channel is marginaly stable 
and then frequently collapses. As a conclusion, Qb — Q* ■ 



P. W. Hammer, N. Piatt, S. M. Hammel, J. F. Heagy, 

and B. D. Lee , Phys. Rev. Lett. 73, 1095 (1994). [21] 

F. Rodelsperger, A. Cenys, and H. Banner, Phys. Rev. 

Lett. 75, 2594 (1995). [22 
D. L. Feng et al, Phys. Rev. E 58, 3678 (1998). [23 
T. John, R. Stannarius, and U. Behn, Phys. Rev. Lett. 
83, 749 (1999). [24 
A. Vella et al, Phys. Rev. E 67, 051704 (2003). 
R. M. Dickson et al. Nature (London) 388, 355 (1997). [25 
X. Brokmann, J.-P. Hermier, et al., Phys. Rev. Lett. 90, 
120601 (2003). [26 

G. Zumofen, J. Hohlbein and C. G. Hubner, Phys. Rev. 

Lett. 93, 260601 (2004). [27 
C.-C. Lo, L. A. Nunes Amaral, S. Havlin, et al. Europhys. 
Lett. 57, 625 (2002). [28 

H. Fujisaka and T. Yamada, Prog. Theor. Phys. 74, 918 
(1985). [29 
N. Piatt, E. A. Spiegel, and C. Tresser, Phys. Rev. Lett. [30 
70, 279 (1993). 

J. F. Heagy, N. Piatt, and S. M. Hammel, Phys. Rev. E [31 
49, 1140 (1994). 

Y. Pomeau and P. Manneville, Commun. Math. Phys. [32 
74, 189 (1980). 

H. G. Schuster, Deterministic Chaos, 2nd Edition (VCH, [33 
New York, 1989). 

F. Zoueshtiagh and A. Merlen, Phys. Rev. E 75, 056313 [34 
(2007). [35 
L. Gostiaux, H. Gayvallet and J.-C. Geminard, Granular [36 
Matter 4, 39 (2002). 

G. Varas, V. Vidal and J.-C. Geminard, Phys. Rev. E 

79, 021301 (2009). [37 
M. Ripepe, S. Ciliberto, M. Delia Schiava, J. Geophys. 
Res. 106, 8713 (2001); M. Ripepe, A.J.L. Harris, R. [38 
Carniel, J. Vole. Geotherm. Res. 118, 285 (2002). 

H. M. Gonnermann and M. Manga, Annu. Rev. Fluid [39 
Mech. 39, 321 (2007). 

S. L. Webb and D.B. Dingwell, J. Geophys. Res. 95, 695 [40 



(1990). 

Y. Lavallee, K.-U. Hess, B. Cordonnier, and B. D. Ding- 
well, Geology, 35, 843-846 (2007). 
E. A. Parfitt, J. Vole. Geotherm. Res. 134, 77 (2004). 
Y. Lavallee, P. G. Meredit, B.D. Dingwell et al. Nature, 
453, 507 (2008). 

L. Caricchi, D. Giordano, L. Burlini, P. Ulmer and C. 
Romano, Chem. Geology, 256, 158 (2008). 
C. Jaupart and S. Vergniolle, Nature 331, 58 (1988); J. 
Fluid. Mech. 203, 347 (1989). 

S. Okumura, M. Nakamura and A. Tsuchiyama, Geo- 
phys. Res. Lett. 33, L20316 (2006). 

S. Okumura, M. Nakamura, A. Tsuchiyama, T. Nakano 
and K. Uesugi, J. Geophys. Res. 113, B07208 (2008). 

A. Tufaile and J. Sartorelli, Phys. Rev. E 66, 056204 
(2002). 

I. L. Kliakhandler, Phys. Fluids 14, 10 (2002). 

H. Pinkerton and G. Norton, J. Vole. Geotherm. Res. 68, 

307-323 (1995). 

T. Divoux, V. Vidal, F. Melo and J.-C. Geminard, Phys. 
Rev. E 77 056310 (2008). 

S. Daugan, L. Talini, B. Herzhaft and C.AlIain, Eur. 
Phys. J. E 7, 73 (2002); Eur. Phys. J. E. 9, 55 (2002). 

B. A. Berg and R.C. Harris, Comput. Phys. Commun. 
179, 443448 (2008). 

E. Bertin and F. Bardou, Am. J. Phys. 76, 630 (2008). 
A.W. Woods and T. Koyaguchi, Nature 370, 641 (1994). 
A. Tzounakos, D. G. Karamanev, A. Margaritis and 
M. A. Bergougnou, Ind. Eng. Chem. Res. 43, 5790 
(2004). 

H.A. Barnes, J. Non-Newtonian Fluid Mech., 70 1-33 
(1997). 

S. T. Thoroddsen, T. G. Etoh and K. Takehara, Phys. 
Fluids 19, 042101 (2007). 

F. H. Zhang and S. T. Thoroddsen, Phys. Fluids 20, 
022104 (2008). 

Please note that the formation of satellite by bubble 



11 



pinch-off and bubble coalescence has been studied in the 
case of a Newtonian fluid and is still an active field of 
research [3^ . [39l IZlj ; This topic still remains to be fully 
characterized in the case of non-Newtonian fluids. 



[41] J. C. Burton and P. Taborek, Phys. Rev. Lett. 101, 

214502 (2008). 
[42] J. Eggers, Rev. Mod. Phys. 69, 865 (1997). 



